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Abstract 

The magnetic catalysis of discrete chiral symmetry breaking in the 2+1 
dimensional Nambu-Jona-Lasinio model is analyzed. A particular atten- 
tion is paid to a possible application of the effect in solid state physics. 
The fermion contribution to the thermal conductivity as a function of the 
dynamical mass (energy gap in the spectrum) is estimated and is shown 
to be suppressed when a non-zero order parameter develops. 



*On leave of absence from Bogolyubov Institute for Theoretical Physics, Kiev 252143, 
Ukraine 



1 



Chiral symmetry breaking induced by an external magnetic field has recently 
been studied in a few 2+1 and 3+1 dimensional relativistic field theories pi 0, |l 
[| H (for some earlier results on the subject see Rcfs. || 0, ||). In particular, a 
magnetic field enhances fermion condensation in the Nambu- Jona-Lasinio (NJL) 
model in 2+1 dimensions. The solution of the gap equation indicates that 
the critical 4-Fermi coupling required to generate a chiral symmetry breaking 
condensate with ("ipip) ^ goes to zero as soon as the magnetic field is switched 
on. In fact, there are two branches of the solution depending on whether the 
effective 4-Fermi coupling g is less than or greater than g c , where g c is the 
critical coupling in the NJL model at zero magnetic field. In the subcritical 
region g < g c the non-trivial symmetry breaking solution disappears as the field 
goes to zero. On the other hand, at g > g c the solution goes to the usual 
symmetry breaking solution in the NJL model without a magnetic field. 

In this letter we promote the hypothesis that this magnetic catalysis of chiral 
symmetry breaking can be used to explain recent experimental results of K. Kr- 
ishana et al |9). They observed a magnetic field induced phase transition in a 
high-T c superconductor (below T c ). The experiment shows that the thermal 
conductivity as a function of temperature and magnetic field displays a sharp 
break in its slope. There are several competing arguments which attempt to 
explain this behavior. In one line of reasoning, the behavior is accounted for 
the formation of an energy gap at the Fermi surface. In d-wave superconduc- 
tors there is no energy gap in the spectrum of the quasi-electrons. The Fermi 
surface consists of four independent degeneracy points. The gap (as well as the 
temperature of the transition) as a function of field strength follows a square 
root law over a large range of fields. We present this as an alternative to the 
work in jlO| which deals with some general scaling properties of a model in a 
magnetic field, the work of Laughlin which relates the phase transition to 
the appearance of an additional parity and time-reversal symmetry violating 
superconducting order parameter which may develop under specific conditions 
[ |12| and other speculations that the behavior is a signal of a phase transition 
in the lattice of magnetic flux lines. While our work was in progress, we be- 
came aware of a paper by K. Farakos and N.E. Mavromatos Jl3) where a similar 
attempt of applying the magnetic catalysis to the explanation of the phase tran- 
sition observed by Krishana et al || was made. They considered the magnetic 
catalysis in a gauge theory with a continuous symmetry, while we study the NJL 
model with discrete symmetry. We also estimate the fermion contribution to 
the thermal conductivity as a function of the dynamical mass, and show that it 
is suppressed when the order parameter develops a nonzero expectation value. 

In our hypothesis, the additional order parameter which signals the phase 
transition is a neutral condensate of quasi- electron-hole pairs around the nodes 
on the Fermi surface in a d-wave superconductor. The distinctive feature is the 
neutrality of the condensate, and the parity as well as time-reversal symme- 
try conserving character of the order parameter (same as in |l^| ). The phase 
transition itself is of second order. 
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The NJL model is a natural candidate for the description of the interactions 
of quasi-electrons in a d-wave superconductor. It is well-known that, when the 
Fermi level is concentrated at points, the spectrum of the quasi-electrons with 
energies close to the Fermi level is linear in the momentum E{k) ~ |fc| and can be 



described by a Dirac Hamiltonian |13L 14 JL5l. The effective action from which 



the Dirac Hamiltonian can be derived is the Lorentz invariant Dirac action with 
two species of 4-component spinors. The two species arise from the existence of 
four inequivalent degeneracy points on the Fermi surface. The Dirac action has 
apparent continuous chiral and flavor symmetry which are not a symmetries of 
the original system which is defined on a lattice. 

The interaction term is a 4-Fermi operator. The choice of suitable operators 
is dictated by the symmetries of the theory. Though the lattice theory does not 
have continuous chiral symmetry, it does have a symmetry which corresponds 
to discrete chiral symmetry — invariance under translations by one site on the 
lattice. Translation by one site interchanges the degeneracy points of the Fermi 
surface and corresponds to a discrete chiral transformation in the continuum. 
For this reason it is natural to choose a 4-Fermi interaction which breaks the 
continuous chiral symmetry to its discrete Z2 subgroup. 

It is also natural to choose a 4-Fermi operator which is invariant under parity, 
time reversal and charge conjugation invariance. Parity and time reversal are 
symmetries of the theory in the absence of an external magnetic field and charge- 
conjugation invariance is equivalent to the approximate particle-hole symmetry 
of the low energy excitations. 

Even with these restrictions, there are both Lorentz invariant and non- 
invariant 4-Fermi operators. The rationale for choosing a Lorentz invariant 
operator is that, it is a combination of Lorentz invariance and chiral symmetry 
which forbids an energy gap in the spectrum of the continuum model. If we 
chose a non Lorentz invariant 4-Fermi operator, radiative corrections would gap 
the quasi-electron spectrum even without an external magnetic field or spon- 
taneous symmetry breaking. A candidate for generating 4-Fermi interactions 
of the kind that we choose, —g (ipi/A is the effective interaction mediated by 
optical phonons. 

Consider the Nambu-Jona-Lasinio model in 2+1 dimensions [l(|. In 
2+1 dimensions the minimal representation of the Dirac algebra is with 2 x 
2 matrices and the minimal Dirac spinor has two components. However, for 
two component fermions, the mass operator ipip is a pseudo-scalar and massive 
fermions violate parity. If there are two species of fermions, the mass operator 
V'i'01 — V'2'02 is a scalar if, besides the spacetime transformation, the fermion 
species are interchanged under parity. The other combination, tpi^i + ip2tp2 is 
a pseudoscalar. It is convenient to describe the two species of fermions using a 
reducible 4x4 representation if the Dirac algebra JItJ . The 7- matrices are: 

03 \ if iai \ 2 _ f 
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where at are the Pauli matrices. 

Thus, the Lagrangian density of our model is given by 



£=i[*,i 7 ^*]+|(^) 2 , i_>) 



where the covariant derivative includes only the vector potential for the external 
magnetic field: 



= 0„ - leA™\ Af = ( 0, -|* 2 , 1^,0 ) . ( :■{ ) 



To apply the 1/N perturbative expansion we assume that fermions have an 
additional flavor index i = 1, 2, . . . , N. In the d-wave superconductor N=2. 

The standard mass term, m^&\&, as is easy to check, is not invariant under 
the discrete chiral symmetry 

* -> 75*, * -> -*7 5 . (4) 

So, it is forbidden in the Lagrangian density. However, due to the magnetic 
catalysis jjj the discrete symmetry is broken dynamically and as a result, fermions 
get a dynamical mass. Let us stress that symmetry breaking appears already 
at an arbitrary small 4-fermion coupling constant. 

The latter fact can be understood by the following argument. Consider the 
Dirac Hamiltonian in a background magnetic field 

I) 

where 

h = ia 2 D 1 - ia 1 D 2 . (6) 

In a background field with magnetic flux $ = J Bd 2 x an index theo- 
rem (an open space version of the Atiyah-Singer Index theorem) and a vanish- 
ing theorem imply that there are [|<E>|] — 1 solutions of the equation h<f> = 
where [|$|] is the largest integer less than |$|, and their chirality is given by 
(T 3 <fi = sign($)([|$|] — 1). Existence of these zero modes of the Hamiltonian 
implies degeneracy of the second quantized ground state - the degenerate states 
correspond to all of the ways that the zero modes could be populated by quasi- 
electrons. Overall charge neutrality implies that half of the zero modes should 
be populated and half of them should be empty. There are chirally symmetric 
populations of the zero modes where half of the modes with four component 



eigenfunctions of H , * + = I ^ j and half of the modes with wavefunctions 
\&_ = ( ? ) are occupied. The chiral condensate is (ty^"/ ^) , and the average 
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condensate (J cPx^^fy) vanishes. There are also chirally asymmetric popula- 
tions of the zero modes, for example, when all of the positive helicity modes X P+ 
arc filled and the negative helicity modes ^_ are empty. They are degenerate 
with the chirally symmetric ones. 

The degeneracy is resolved by adding the interaction terms in the Hamilto- 
nian and diagonalizing the matrix of first order corrections. In this case, the 
perturbation given by the 4-Fermi operator — G/2 ('i"!') resolves the degener- 
acy by lowering the energy of the states with maximally chirally asymmetric 
populations. It is also easy to see that a similar argument would not work when 
there is no magnetic field since in that case there are no zero modes and the 
density of states at the Fermi surface vanishes like \E — Ep\. The excited states 
which would have to be populated in order to have a non-vanishing condensate 
have positive energy and therefore there is a finite critical coupling which g must 
achieve before the 4-Fermi can lower the energy of the chirally asymmetric state 
below that of the ground state. 

Conventionally, the symmetry breaking is found using the effective potential 
in the large N expansion. The analysis of the model with the Lagrangian density 
in (H) repeats, up to minor differences appearing due to the absence of Nambu- 
Goldstone bosons, the analysis in pL We shall skip all of the intermediate 
calculations and go directly to the gap equation 



v^l + ^H + Oa/A). (7) 



iTfldyn 2 2 



Here we introduced the mass scale parameter [i = 2A(g c — g)/gg c (by definition, 
g = NGA/ir and g c — ^Jtt). In addition, the magnetic length I = 1/ ^J\eB\ is 
used instead of the magnetic field strength B. From now on we assume that the 
value of cutoff A is much bigger than all other parameters in the model, so that 
further we can neglect all corrections of order 1/A or higher. 

The gap equation (Q) is written for the "subcritical" value of coupling con- 
stant, g < g c (where g c is the critical value of coupling in the model without 
the magnetic field). If we try to go into the supercritical region the parame- 
ter changes the sign. In this paper, however, we are not going to study the 
case when g > g c . As is clear, by fixing the dimensionful parameter /z, one 
already defines the NJL model. As soon as the fj, is fixed, we can talk about 
strong or weak external magnetic fields applied to the system. The quantitative 
characteristics of that is the value of dimensionless parameter fil. Indeed, the 
vanishing field will correspond to infinitely large value of fxl. While increasing 
the field, the parameter /il gets smaller, and eventually goes to zero when field 
becomes infinitely strong. Of course, different values of ill can also be inter- 
preted from another viewpoint: small parameter id corresponds to near-critical 
value of coupling constant, while large value of ill corresponds to weak coupling. 

For our analysis below, the first approach seems to be more appropriate since 
we are interested in the dependence of dynamical mass on the magnetic field 
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strength without changing the "parent" NJL model. It would also allow us to 
make the analysis almost independent of a specific value of coupling constant. 

So, let us consider, two limiting cases of field strength: (i) vanishingly weak 
field (or rather y/\eB\ <C zx), and (ii) infinitely large field (or just ^J\eB~\ ^> /i). 

To get the gap equation in the case of weak field, one has to notice that 
large positive value in the right hand side of (0) is obtained for small values of 
fndynl when the first term dominates over the zeta function term. Thus, the 
approximate gap equation reads 

^ — -j — - — i (8) 

so that the solution for the dynamical mass is proportional to the magnetic field 
strength, 

\eB\ 

Ttldyn = ■ (9) 

M 

Here we would like to note that a simple numerical analysis of the gap equation 
(0) shows that the solution is indeed close to this asymptotic behavior up to a 
few percent but only in the region where the ratio y/\eB\/fx is less than around 
0.005. If the ratio gets bigger than that deviations are more prominent. 

The case of strong magnetic field is slightly more complicated. The first 
approximation for the gap equation in this case is obtained by substituting zero 
in the left hand side of the gap equation (^) : 



V2C(il + ^). (10) 



Iwidyn ^ ^ 

This equation can easily be solved numerically, leading to the result rridynl = 
cq Ri 0.446. And the solution for the dynamical mass is proportional to the 
square root of the field, 



rridyn — cq VjeSf + ciAt, (11) 

where the second term (with c% = —0.173) is the next to leading term in Taylor 
expansion around [il = 0. It appears when the magnetic field strength is large 
but finite. Even though the solution in (nil) was obtained under the assumption 



\/\eB\ ^> fi, it should be rather reliable even for values of ^/|e_B| of order //. In 
fact, the numerical solution of the gap equation, given in Figure 1, confirms our 
statement. More than that, as a simple analysis shows, the numerical solution 
can be well approximated by a the square root dependence in the whole region 
of magnetic fields except for a small region of weak fields. 

Thus, summarizing, the dependence of the dynamical mass on the exter- 
nal magnetic field is linear for very small fields (small in comparison with the 
inherent scale /x), then the dependence gradually approaches the asymptote in 



G 



Eq. ( |l 1| ) when the square root of field gets bigger than fi. The graphical solution 
is shown in Figure 1. 

In view of application the result to the phase transition mentioned in the 
introduction, we have to require that the actual value of the parameter /i is 
small enough in comparison with the square root of the magnetic field while 
given in the appropriate units. So, that the small linear region of the dynamical 
mass dependence is not experimentally seen. This assumption is not unnatural 
at all if we recall that the asymptotic behavior in ( |TT| ) has a rather wide range 
of validity, while the linear dependence appears only for really very weak fields. 

The gap equation at a finite temperature can also be written right away 
|Q. However, in that case it becomes really complicated and we are not going 
to write it down here. The qualitative picture though is clear enough even 
without solving the equation. In particular, we know that the magnetic field 
helps symmetry breaking, while temperature works against it. As a result, we 
have a competition of those two influences. If we keep the value of the field 
being constant, the critical temperature at which the order parameter (gap 
in the energy spectrum) disappears is, up to an unimportant numerical factor, 
proportional to the value of the dynamical mass at zero temperature md yri {\eB\) 
[0- H- [01 ■ On the other hand, considering the mass as a function of the field, we 
can invert the equation for the critical temperature, and obtain the critical value 
of the magnetic field at a given temperature. So, looking for the dependence of 
the dynamical mass on the field strength at a finite temperature, we will find 
that the gap stays zero until the value of the field is less than critical, and after 
field becomes stronger than critical, a non-zero gap develops and grows with 
increasing the field. Later we will make use of this qualitative picture. 

We would like to mention that one may be tempted and try to use the 
magnetic catalysis of 2+1 dimensional quantum electrodynamics instead of the 
NJL model. It is indeed a very attractive idea since, in a rather physical limit 
of a weak magnetic field (and apparently in range of validity of the large N 
expansion), the dynamical mass as a function of the magnetic field follows the 
square root law dependence [EOl : 



m dyn ~ ln(4irv N), (12) 

where v$ w 0.14. Note that this dependence is quite natural from the following 
point of view. As is well known, in 2+1 dimensional QED it is 1/N rather than 
e 2 that plays the role of the coupling constant (lj, |||] . The value of e 2 is more 
like an ultraviolet cutoff in the model. Therefore, one effectively is left with a 
single infrared dimcnsionful parameter, namely -\/]eB\, and the dynamical mass 
should be proportional to this only one dimcnsionful parameter. 

While considering quantum electrodynamics one should remember, however, 
that breaking the continuous symmetry is forbidden in 2+1 dimensions at any 
finite temperature. In order to avoid this difficulty one should apply additional 
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assumptions, like those in |L3| where the realization of the BKT phase is con- 
jectured. 

Another thing that can obscure the effect of magnetic catalysis in 2+1 di- 
mensional QED is the parity violating mass term that seems to be generated 
due to the induced Chern-Simons photon mass in a magnetic field p3fl . Even 
if the effect of magnetic catalysis does not disappear, the quantitative analysis 
may somewhat change. 

It is possible to give a simple estimate of the thermal conductivity depen- 
dence on the magnetic field in the model at hand. 

To obtain the expression for the thermal conductivity we apply the familiar 
method of linear response ^4|. There are, however, a few things that we need 
to mention. First, we work with zero chemical potential, so that the general 
expression for the conductivity is given in terms of the energy-current correlation 
function alone. And second, to get a finite thermal conductivity we will need 
to modify the free fermion Green function in the magnetic field by introducing 
an effective lifetime of quasiparticles. The latter can be interpreted as a result 
of fermion scattering at impurities of a superconductor, which is always present 
in real samples. Such an approximation is presumably good enough for our 
purposes, since the only thing we are interested in is a change in the value of 
the thermal conductivity that results from the magnetic catalysis of symmetry 
breaking and generation of the energy gap. 

As we mentioned above, the expression for the thermal conductivity can be 
given in terms of the energy-current (momentum) correlation function pij 

oo 

<M = wl dt l dXTr {<°° pi (°) p ^ + lA )l e ~ ia "< ( 13 ) 



where V is the volume of the system, 

P l (0) = % - J d 2 x (^ 7 °<9V " 0V7V) , ( 14 ) 
pi(t) = e tHt P* \0)e- lHt (15) 
and po is the equilibrium density matrix: 

P0 = V^, Z = Tr(e-^). (16) 

Notice here that the action is invariant under translations even in the presence 
of a constant external magnetic field, and that the definition of the momentum 
coincides with that for the theory without an external field. Here we do not 
include the contribution of the gauge fields to the momentum since we are 
interested in the fermion component of the thermal conductivity. 
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By using the definition in Eq. (|l3|), one can obtain the following expression 
for the static thermal conductivity of an isotropic system pjl |2(1 [27|] 



K Cl - 

~ TV 



OO 

l -Im J tdtTr{p P\t)P\Q)}. (17) 



Let us introduce the following thermal Green function: 



G(t) = ±Tr(p e TH P*e- TH P l ), (18) 



G{iu m ) = f G(r)e w ^ T dr, (19) 



where v m = 2-KTm. It is this Green function G(v m ) that can perturbatively be 
calculated in a theory. The simplest Feynman diagram gives 

G(p = 0, iu m ) =iT 7^ khr (l° S ( iLJ n, k)ry S(iu n + iv m , k)) , (20) 

It turns out that the thermal conductivity can be given by the discontinuity of 
this Green function (see (2^. 26, 27 for details) 

K el = — lim - (G(p = 0,iv m = uj + i0+) - G(p = Q,iv m = u)-iQ + )) . (21) 
AT ui^a lu v ' 

To calculate the Green function in Eq. ( p0[ ) , it is convenient to use the following 
spectral representation for the fermion thermal Green function j2(| [27| 

+ 00 

S(i. n ,k)= [ £^1*1, ( 22 ) 

J Z7T lLU n — LU 

— oo 

where a(u),k) = 2ImS{iu) n = lu — i0 + ,k). After substituting this spectral 
representation into Eq. (p0|), one can easily perform the sum over the Matsubara 
frequencies. Then the result reads 



i 

G(p = 0, iv m ) = 



d 2 k [ duj-L f duj 2 k i smh('^=^ 



IT 



2 J (2tt) 2 7 2vr J 2ir cosh cosh 
tr ( 7 a(^i,fc)7 a(^ 2 ,fc)) 
(lox -u> 2 - iv m ) 



(23) 



As we see, this representation is indeed very convenient for extracting the discon- 
tinuity that is needed for calculation of the thermal conductivity [see Eq. (pT 



First, let us calculate the thermal conductivity in the lowest Landau level ap- 
proximation without taking into account any interactions in the theory. The 
corresponding fermion Green function is given by 

s ^- k) = exp (rwu (1 - 11 1 ] ' (25) 

so that the spectral density reads 

a(u), k) = 27rsign(w) exp (- 5 (uj 2 - to 2 ) [w 7 ° + to] (l - ij 1 -/ 2 ) . (26) 

Making use of this explicit expression for the spectral density, one obtains the 
thermal conductivity right away 



el 



2 



K 



\eB 



16T 2 V2T 



sech 2 — U(0). (27) 



This result looks meaningless because of the (S-function in the right hand side 
which is infinite. Note, however, that this infinity appears due to the 5-like peak 
of the spectral density a(k, uj) in Eq. (|26|). The latter is a direct consequence of 
our dealing with a free theory. In practice, the spectral density of any physical 
interacting model is always spread over a finite region of energies. 

Instead of actual consideration of the interaction effects on the spectral den- 
sity in our model, we just modify the latter by applying simple phcnomenological 
arguments. In particular, we introduce the width of quasiparticle states (or in 
different words, the inverse quantity of the lifetime), r = 1/r <C to, according 
to the following rule 



S (uj 2 — to 2 ) = - — (5 (uj + to) + S (uj — to)) 
2to 

i / r 



2m7r \(uj + to) 2 + T 2 (uj- to) 2 + r 2 



(28) 



At this point we may not specify the nature of the quasiparticle lifetime, but use 
it as a phenomenological parameter. However, we note that such broadening of 
the spectral density can be accounted for the scattering of fermions off impurities 
which always present in real systems [^6[ ^71 . 

Now we are in a position to estimate the thermal conductivity in our modified 
theory. In order to preserve the validity of our symmetry breaking analysis due 
to the magnetic catalysis, we assume that the width T is considerably smaller 
than the dynamical mass, T <C to. Then substitution of the spectral density 
into the expression for the conductivity ( p4| ) leads to the final result 

„/ \eB\ 2 . o / to 



sech 2 ^). (29) 



327rT 2 r V277 
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Let us apply this simplified result to the study of thermal conductivity depen- 
dence on the magnetic field at a fixed temperature. Qualitatively, we get the 
following picture. While the magnetic field is weaker than the critical value at a 
given temperature, there is no energy gap in the spectrum of the model. So, we 



have to substitute m(T, \ eB\ < \eB c \) = in Eq. (27). As a result, we observe 



that the thermal conductivity increases with the magnetic field: 

\eB\ 2 

Kel = 32^T¥' (30) 

This seems to be just the opposite what the experiment sees However, 
we believe that this increase might be the artifact of inappropriate use of the 
approximation (see below) as well as a result of absence of the Fermi surface in 
our model and should not be taken as a real prediction. 

On the other hand, after the magnetic field gets bigger than the critical 
value, \eB\ > \eB c \, the energy gap m(T, \eB\ > \eB c \) ^ develops and grows 
with increasing the field strength. So, that gradually, the dependence of the 
thermal conductivity becomes exponentially damped by the energy gap: 

We expect that this damping of the conductivity is a rather general feature 
which should be present as soon as an energy gap is developed in the model. 
We remind that something similar also happens with the thermal conductivity 
when an ordinary superconducting order parameter develops [ p6| . 

In support of the estimates made above, let us comment on the lowest Lan- 
dau level approximation that was used in the derivation of Eq. ( p7j ) . Obviously, 
this approximation could be reliable only if we stay in the region where the dy- 
namical mass is much smaller than the scale of the magnetic field. Let us remind 
that this condition was indeed satisfied in the case of the magnetic catalysis in 
QED at zero temperature. While considering finite temperatures, we have to 
require also that the temperature is much less than the square root of the field 
strength. This latter requirement gives us one more argument in support of 
reliability of result for the thermal conductivity (|l]) in the case of supercritical 
fields, while adds more distrust in the analysis for subcritical magnetic fields. 

Concerning the NJL model, we mention that the ratio m/\/\eB\, as is seen 
from Eq. (|TT1), approaches 0.446 in the strong field limit. It is this region that 
reproduces the square root dependence of the thermal conductivity as a function 
of the magnetic field. Consider, however, that approximately the same depen- 
dence law remains valid down to the values of the ratio m/y/\eB\ as small as 
0.1. Therefore, even though the lowest Landau level approximation may not be 
perfect in case of the NJL model, it should be reliable enough for the thermal 
conductivity estimate and for the conclusion about the damping of the latter in 
the supercritical region. 
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The study of the 2+1 dimensional relativistic NJL model, modeling the 
dynamics of quasiparticles close to the nodes of the Fermi surface in d-wave 
superconductors, reveals the dependence for the mass as a function of the mag- 
netic field very close to square root law for rather general assumptions. The 
appearance of the energy gap is viewed here as a result of breaking a discrete 
symmetry due to neutral condensate of fermion-antifcrmion pairs. In condense 
matter language such a situation could be interpreted as the condensation of 
electron- hole-like quasiparticle pairs. The results of our paper give an alterna- 
tive explanation to the recently observed |J phase transition in cuprites below 
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